Declaration 



Honorable Commissioner of Patents: 
I, Xu Youhao, hereby declare that: 

1. I'm living in Beijing, People's Republic of China. 

2. I'm a citizen of People's Republic of China. 

3. From September 1981 to July 1985, I studied in the Department of 
Petroleum Refining in East China University of Science and Technology and 
obtained my Bachelor's degree; from September 1985 to July 1987, 1 studied in 
the Department of Chemical Engineering in East China University of Science 
and Technology and obtained my Master's degree; from September 2001 to 
December 2006, 1 studied in the Research Institute of Petroleum Processing and 
obtained my Doctor's degree. 

4. Since August 1987, I've been devoting myself to the development of 
catalytic cracking technology in the Research Institute of Petroleum Processing, 
and have twenty years' experience. I took part in the development of Deep 
Catalytic Cracking (DCC) technology, and I'm currently in charge of the 
development and application of a newly series of FCC processes for producing 
maximum isoparaffms in cracked naphtha or clean gasoline and propylene (MIP 
and MIP-CGP). I have written more than forty papers on this subject, six of 
which were published in internationally significant periodicals and international 
academic conferences, and applied for, altogether, more than seventy patents 
domestically and abroad. 

5. In US Patent Application 09/553,990, I propose a novel riser reactor, 
wherein the ratio of the diameter of a second reaction zone to the diameter of a 
first reaction zone ranges from 1.5: 1 to 5: 1. However, in catalytic cracking of 
the prior art, the ratio of the diameter of a second reaction zone to the diameter 
of a first reaction zone in a riser reactor should be less than 1.5: 1. 

6. On the basis of the European patent EP0171460, Kmecak et al., the 
Examiner in charge of this application drew the conclusion, after manually 



making a measurement in Figure 8 of EPO 171460, that "DR 2 : DRi is equal to 
approximately 3: 1". In my point of view, this is impossible because of the 
following reasons: 

6.1. In catalytic cracking of the prior art, the ratio of the diameter of a 
so-called second reaction zone to the diameter of a so-called first reaction zone 
in a riser reactor should be less than 1.5: 1. 

6.2. Figure 8 in EP0171460 is merely a schematic drawing without the 
specific dimensions available. Thus, the value obtained by manual measurement 
on the basis of the schematic drawing is inaccurate. One skilled in the art would 
never do this. 

6.3. If "DR 2 : DRi is equal to approximately 3: 1", it will be obtained upon 
scientific calculation (see Items 8 and 9 for details) that the ratio of the velocity 
at the outlet of the first reaction zone to the velocity at the inlet of the first 
reaction zone, i.e. Ui/u 0 , is 18.7, which is far greater than the generally 
recognized ratio of less than 3 in the art. Thus, the Examiner's conclusion 
cannot be fulfilled. 

7. In my point of view, the present invention achieves unexpected technical 
effects, which can be seen clearly from Example 1 , Example 2, Example 3 and 
the corresponding Comparative Example, as given herein below. 

The Comparative Example A and Example 3 are complementary tests. 
Example 1 and Example 2 are recited in the original application document. The 
properties of feedstock used in the Examples and the Comparative Example A 
and the properties of the catalyst are respectively listed in Tables 1 and 2. 

Example 1 

This Example showed the use of catalyst A in a novel pilot plant riser 
reactor according to the present invention for the production of isobutane and 
isoparaffm-enriched gasoline by means of catalytic conversion in the riser 
reactor. 

The total height of the prelift zone, the first reaction zone, the second 
reaction zone, the outlet zone of the riser reactor was 15 meters, in which the 



height of the prelift zone with the diameter of 0.025 meter was 1.5 meters, the 
height of the first reaction zone with a diameter of 0.025 meter was 4 meters, 
the height of the second reaction zone with a diameter of 0.1 meter was 6.5 
meters, the height of the outlet zone with a diameter of 0.025 meter was 3 
meters. The isotrapezia vertex angle of the vertical section of the conjunct 
section between the first reaction zone and the second reaction zone was 45°. 
The isotrapezia base angle of the vertical section of the conjunct section 
between the second reaction zone and the outlet zone was 60°. 

The operating conditions and product slate are listed in Table 3, and the 
gasoline properties are listed in Table 4. 

Example 2 

This Example showed the use of gasoline with high olefin content as a 
quenching medium in a novel riser reactor for the production of isoparaffm 
enriched gasoline by means of catalytic conversion in the riser reactor. 

The total height of the prelift zone, the first reaction zone, the second 
reaction zone, the outlet zone of the riser reactor was 15 meters, in which the 
height of the prelift zone with the diameter of 0.025 meter was 1.5 meters, the 
height of the first reaction zone with a diameter of 0.025 meter was 4 meters, 
the height of the second reaction zone with a diameter of 0.05 meter was 6.5 
meters, the height of the outlet zone with a diameter of 0.025 meter was 3 
meters. The isotrapezia vertex angle of the vertical section of the conjunct 
section between the first reaction zone and the second reaction zone was 45°. 
The isotrapezia base angle of the vertical section of the conjunct section 
between the second reaction zone and the outlet zone was 60°. 

The catalyst and gasoline feedstock used in this Example were the same as 
those used in Example 1 . The gasoline produced in comparative example 1 of 
the original description as a quenching medium was charged into the conjunct 
section between the first reaction zone and the second reaction zone. The other 
steps were basically the same as those in Example 1 . 

The operating conditions and product slate are listed in Table 3, and the 
gasoline properties are listed in Table 4. 



Example 3 



This Example showed the use of gasoline with high olefin content as a 
quenching medium in a novel pilot plant riser reactor according to the present 
invention for the production of isobutane and isoparaffin enriched gasoline by 
means of catalytic conversion in the riser reactor. 

The total height of the prelift zone, the first reaction zone, the second 
reaction zone, the outlet zone of the riser reactor was 15 meters, in which the 
height of the prelift zone with the diameter of 0.025 meter was 1.5 meters, the 
height of the first reaction zone with a diameter of 0.025 meter was 4 meters, 
the height of the second reaction zone with a diameter of 0.0375 meter was 6.5 
meters, the height of the outlet zone with a diameter of 0.025 meter was 3 
meters. The isotrapezia vertex angle of the vertical section of the conjunct 
section between the first reaction zone and the second reaction zone was 45°. 
The isotrapezia base angle of the vertical section of the conjunct section 
between the second reaction zone and the outlet zone was 60°. 

The catalyst and gasoline feedstock used in this Example were the same as 
those used in Example 1 . The gasoline produced in comparative example 1 of 
the original description as a quenching medium was charged into the conjunct 
section between the first reaction zone and the second reaction zone. The other 
steps were basically the same as those in Example 1 . 

The operating conditions and product slate are listed in Table 3, and the 
gasoline properties are listed in Table 4. 

Comparative Example A 

Compared with the previous Example 2, this comparative example was 
conducted under the same condition except that the diameter of the second 
reaction zone is 0.03 meters and thus the diameter ratio of the second reaction 
zone to the first reaction zone is 1 .2: 1 . 

The operating conditions and product slate are listed in Table 3, and the 
gasoline properties are listed in Table 4. 

It can be seen clearly from Tables 3 and 4 that: 



1) In Example 1, 5.65wt% of LPG was isobutane; 
In Example 2, 5.98wt% of LPG was isobutane; 
In Example 3, 5.62wt% of LPG was isobutane; 

however In Comparative Example, 3.02 wt% of LPG was isobutane (see 
Table 3); 

2) In Example 1, Iso-paraffins in gasoline was 36.00wt%; 
In Example 2, Iso-paraffins in gasoline was 43.86wt%; 
In Example 3, Iso-paraffins in gasoline was 34.8wt%; 

however In Comparative Example, Iso-paraffins in gasoline was 
21.33wt%(seeTable 4); 

In other words, when the ratios of the diameter of the second reaction zone 
to the diameter of the first reaction zone are 4:1, 2:1 or 1.5:1 respectively, the 
present invention has a significant advantage in respect of effects over the prior 
art in which the corresponding diameter ratio is 1 .2: 1 . 



Table 1 



Density (20 °C), kg/m 3 


890.5 


Kinematic Viscosity (100 °C), mm 2 /s 


5.08 


Carbon Residue, wt% 


0.7 


Pour Point, °C 


40 


Total Nitrogen, wt% 


0.16 


Sulfur, wt% 


0.53 


Carbon, wt% 


85.00 


Hydrogen, wt% 


12.62 


Heavy Metal Content, ppm 




Ni 


0.16 


V 


0.15 


Fe 


— 


Cu 


— 


Na 


0.45 


Distillation, °C 




IBP 


278 


10% 


385 


30% 


442 


50% 


499 


70% 




90% 




EP 





Table 2 



Catalyst Name 


A 


Trade Mark 


ZCM-7 


Chemical Composition, wt% 




Aluminum Oxide 


46.4 


Sodium Oxide 


0.22 


Ferric Oxide 


0.32 


Arvnarpnt RnlV T~)pnQitv Vcr/m^ 


690 


"Pnrf* \/olnmp ml/a 
-TVJ1C VU1U.111C, 1111/ ii 


0 38 


Surface Area, m 2 /g 


164 


Attrition Index, wt% hr" 1 




Particle Size Distribution, 
wt% 




0-40 microns 


4.8 


40-80 microns 


47.9 


>80 microns 


47.3 



Table 3 





Example 1 


Example 2 


Example 3 


Comparative 
Example A 


the diameter ratio of 
the second reaction 

zone to the first 
second reaction zone 


4:1 


2:1 


1.5:1 


1.2:1 


Operation condition: 




Reaction Temperature, 
°C 


The First Reaction Zone 


545 


545 


545 


545 


The Second Reaction 
Zone 


495 


495 


495 


495 


Reaction Time, sec 


5.0 


5.3 


5.0 


5.0 


The First Reaction Zone 


1.0 


0.8 


1.5 


2.2 


The Second Reaction 
Zone 


3.5 


3.9 


3.0 


2.3 


The Outlet Zone 


0.5 


0.6 


0.5 


0.5 


C/O Ratio 


4.5 


5.0 


4.5 


4.5 


S/O Ratio 


0.05 


0.05 


0.05 


0.05 


Product Slate, wt% 




Dry Gas 


1.83 


1.78 


1.95 


4.65 


LPG 


16.11 


17.51 


16.54 


13.17 


in which Isobutane 


5.65 


5.98 


5.62 


3.02 


Gasoline 


46.86 


47.98 


46.33 


42.59 


LCO 


23.44 


22.30 


23.50 


23.12 


HCO 


7.77 


6.22 


7.65 


10.17 


Coke 


3.88 


4.00 


4.00 


6.00 


Loss 


0.11 


0.21 


0.03 


0.30 



Table 4 





Example 1 


Example 2 


Example 3 


Comparative 
Example A 


the diameter ratio of 
the second reaction 
zone to the first second 
reaction zone 


4:1 


2:1 


1.5:1 


1.2:1 


Density (20 °C), 
kg/m 3 


743.6 


745.3 


742.6 


748.2 


Octane Number 




RON 


90.0 


90.1 


90 


90.5 


MON 


79.0 


80.9 


79 


79.3 


Induction Period, 


>1000 


800.0 


>1000 


343 


J_fA.li>lClll VJU.1I1 9 

mg/ 100ml 


2.0 


2.0 


2.0 


2.0 


Sulfur, wt% 


0.0095 


0.01 


0.009 


0.016 


Nitrogen wt/'o 


0.0028 


0.003 


0.003 


0.004 


WL/O 


86.14 


86.51 


86.15 


86.72 


Hydrogen, wt% 


13.72 


13.42 


13.73 


13.19 


uisiuianon, 




IBP 


46 


48 


46 


48 


10 % 


73 


75 


73 


75 


30% 


95 


97 


95 


96 




114 


118 


114 


120 




143 


144 


143 


145 


90 % 


171 


173 


171 


173 


EP 


202 


203 


202 


200 


Gasoline 

composition, wt%> 




Paraffins 


41.01 


47.87 


40.01 


25.61 


in which Iso-paraffins 


36.00 


43.86 


34.8 


21.33 


Naphthenes 


7.20 


7.45 


7.0 


6.5 


Olefins 


28.11 


20.51 


28.56 


40.69 


Aromatics 


23.68 


24.17 


24.43 


27.20 



8. Directed at the part of "Response to Arguments" on pages 10-16 in the 
Office Action on December 13, 2007, 1 deem 

(1). The density in Fig. 4 of US3246960 is fundamentally different from the 
density in my previous Declaration. The density mentioned in my previous 
Declaration is the density of gas in the riser reactor and the objects mentioned 
during the calculation are all gas, while the density in Fig. 4 of US3246960 is 
the density of mixture of gas and solid in the riser reactor (see column 4, line 
74; column 2, lines 9-11). There are three reasons listed one by one below. 

The first reason : 

column 4, lines 73-74 of US3246960 clearly states that "the density of a 
fluidized mixture", and column 2, lines 9-10 clearly states that " - -the 
suspension has an average density - 

The second reason : 

According to the definition of the mixture of gas and solid (Powder 
Technology 82, 1995, Page 271) 

Pm = Pg X £ + Pp X (l ~ £ ) 

wherin p m represents the density of the mixture of gas and solid, p g represents 
the density of gas, p p represents the density of solid catalyst, s represents 
voidage. Further conversion of the above equation produces : 

P m =P P -(Pp-P g ) X£ 

Because p g «p p , the equation above can be simplified as: 

P m ~P P x i l - £ ) 

The density of catalytic particles is usually 1500 kg/m 3 , the range of the 
voidage in the riser reactor is usually 0.80-0.98 (the above equations and 
data are all common knowledge of a person ordinarily skilled in the art) , so 
the density p m of the mixture of gas and solid is 30-300 kg/m 3 , i.e. 
1.87-18.73 lb./cu.ft., which is consistent with 3.6-27.3 lb./cu.ft. in the right 
table of Figure 4, which proves that the data in Figure 4 of US3246960 is the 
density of the mixture of gas and solid rather than the density of the gas 



mentioned in the present invention. 



The third reason : 

Figure 3-3-11 on pages 128-130, Fluidized Catalytic Cracking Technology 
(by Lu Chunxi and Wang Zhu' an, see Annex Al and its English translation 
Annex A2 of the present Declaration) tests the average density of the mixture 
of gas and solid in the riser under various operating conditions is 46.3-53 
kg/m 3 , i.e 2.89-3.31 lb./cu.ft., which is consistent with that of 1-6 lb./cu.ft. 
in column 2, line 1 1 of US3246960. Said data further proves that the data in 
Figure 4 is the density of the mixture of gas and solid rather than the density 
of the gas mentioned in the present invention. 

Furthermore it can also be seen that the density of the gas (i.e. p! mentioned in 
my previous Declaration) is far less than the density in Figure 4 of US3246960 
(3.6-27.3 lb./cu.ft) from the below calculation (2).l) (i.e. 0.1561b./cu.ft. 
«3. 6-27.3 lb./cu.ft.) . 

(2). The Examiner questions the correctness of p l = p 2 in my previous Declaration. 
Here I evidences based on the data of Chemical Engineering Science 62 (2007), 
Page 1 1 84— 1 198(see ANNEX B) that p, = p 2 is correct. Based on the Figure 4 
on Page 1193, cross section 1 and cross section 2 are respectively taken to 
simulate the two cross sections of First Reaction Zone Outlet and Second 
Reaction Zone Inlet, 




Riser Height (m) Riser Height (m) 

in which CS l=cross section 1 and CS 2=cross section 2. 
1) First Reaction Zone Outlet(cross section 1) 

At the First Reaction Zone Outlet, Pressure P=185800Pa, Temperature 

T=532 °C , and the gas is composed of: 



Gas Components 


Mass Fraction, wt% 


Molar Mass, kg/kmol 


Steam 


5.89 


18 


Fuel Gas (FG) 


2.79 


20 


Liquefied Petroleum Gas (LPG) 


7.86 


50 


Gasoline (GLN) 


33.27 


100 


LCO 


18.41 


200 


Feedstock 


31.78 


500 



The molar mass of the gas at this position can be calculated as 
100 



5.89/18 + 2.79/20 + 7.86/50 + 33.27/100 + 18.41/200 + 31.78/500 
and the density of the gas is 



= 89.9 kg/kmol 



185800x89.9 



= 2.496 kg/m 3 =0.156 lb./cu. ft. 



8.314x(273 + 532)xl000 

2) Second Reaction Zone Inlet ( cross section 2 ) 

At the Second Reaction Zone Inlet, Pressure P=184714Pa , 
Temperature T=530°C , and the gas is composed of: 



Gas Components 


Mass Fraction, wt% 


Molar Mass, kg/kmol 


Steam 


5.89 


18 


Dry Gas 


3.23 


20 


Liquefied Gas 


8.77 


50 


Gasoline 


34.91 


100 


Diesel Fuel 


19.64 


200 


Feedstock Oil 


27.56 


500 



The molar mass of the gas at this position can be calculated as 
100 



5.89/18 + 3.23/20 + 8.77/50 + 34.91/100 + 19.64/200 + 27.56/500 ' 
and the density of the gas is 

= 2.37 kg/m 3 = 0.148 lb./cu. ft. 



85.7 kg/kmol 



184714x85.7 



?.314x(273 + 530)xl000 



It is clear that the density of gas at cross section 1 (pi) is 0.156 
lb./cu.ft which is approximately equivalent to the density of gas at 
cross section 2 (p2) of 0.148 lb./cu.ft. Thus it can be deemed that 
pl=p2. 

(3). Page 12, Lines 10-14 of this Office Action states that " Furthermore, 

the disclosure of Kmecak, in combination with the teachings of Williams, would 
have suggested to one having ordinary skill in the art that the dimensions of the 
zones in the riser could be varied in order to optimize the reaction conditions, 
e.g., residence time, within a particular zone, for establishing the conversion of 
a particular hydrocarbon containing feed stream to the desired products." 

I can not agree with the examiner. I deem that the adjustments to the 
diameter and length of riser (EP 0171460, Kmecak) are very limited other 



than unlimited. 

According to the requirement of production, the structural dimensions of 
riser reactor (including diameter and length) are adjustable in design so as to 
obtain the best product distribution. But since catalytic cracking adopts catalysts 
containing zeolite molecular sieve, and the reactions are very quick, the riser 
reactors used in the current industry imposes stringent requirements on the 
reaction time, which is usually 2-3 seconds. Therefore, the adjustments to the 
diameter and length of riser are very limited other than unlimited. For example, 

1 ) Regarding the riser with the length of 49m and the diameter of 1.2m, if 

the reaction time is 2 seconds, the volume flow is: 

49 1 ? 2 

— x^rx— = 27.7088 m 3 / s 



If the reaction time is changed to 3 seconds, then when the volume flow and 
the length of the riser are the same, the diameter of the riser is: 



Thus, it can be seen that when other conditions do not change, and the 
reaction time is extended from 2 seconds to 3 seconds, the diameter of the 
riser is increased from 1.2m to 1.47m, the range of increase being very 
limited and the ratio of the two being merely 1.225:1. 

2) If the diameter of the riser does not change (1.2m), then when the 
reaction time is extended from 2 seconds to 3 seconds, under the condition 
of the same processing quantity, the height of the riser will increase from 
49m to 73.5m, which is usually not adopted in engineering due to very high 
investment and operating cost. 

It can be seen from the above example that conventional FCC riser 
reactors are restricted by reaction time and thus the range of its diameter 
change is very limited in engineering designs, the expansion ratio of the 
diameters usually being less than 1.5:1. The present invention starts with 
catalytic cracking reaction chemistry and breaks through the intrinsic 
concept mainly based on cracking in conventional FCC reactions. Cracking 
reactions and conversion reactions are conducted respectively in different 
reaction zones, thereby extremely extending the reaction time of catalytic 



2 



4 




cracking reactions up to 3-25 seconds, such that the expansion ratio of the 
diameters in different reaction zones can reach 1.5 : 1-5.0:1. 

(4). Here I detailedly prove that the diameter ratio of the second reaction zone to 
the first reaction zone does not exceed 1.5:1, to say nothing of 3.0:1 in the 
reference EP 0171460(Kmecak) as cited by the Examiner by referring to Fig 1. 



1) . The symbols are defined as follows: 
u — velocity, m/s; 

h — height of the riser, m 

t — residence time, s 

V — volume flow rate, m 3 /s; 

A — cross-sectional area of the riser, m 2 ; 

d — diameter of the riser, m; 

x — ratio of the diameter of the second reaction zone to the diameter of the first 

reaction zone, dimensionless; 

P — pressure, Pa; 

M — molar mass, g/mol; 

R— gas constant, 8.314 J/mol»K; 

T — temperature, K 

p — density, kg/m 3 . 

2) . Symbols t, p and V are defined by the following equations: 

a) 

u 

PXM 

p = r7? (2) 

/d\ 2 

V = U X A = U X TT X ^-J (3) 

3) . Prerequisites or Assumptions for said Deduction 

I. No chemical reaction occurs in the prelift zone, the fluid is in a uniform 
motion state, and the average velocity is represented by u 0 . 

II. The velocity in the first reaction zone exhibits a logarithmic increase 
tendency along the axial direction of the riser, and the average velocity in the 
first reaction zone is represented by the logarithmic mean value of the velocity 
at the inlet (u 0 ) and the velocity at the outlet (ui), which is calculated as 
described hereinbelow: 

Ui-Uq 
lnu-L - lnu 0 



III. The basic reactions in the second reaction zone are conversion 



reactions, and molecular expansion may be neglected. The state of the fluid 
may approximately be regarded as uniform motion state. The average velocity 
is presented by u 2 . 

IV. The velocity in the conjunct section between the first and second 
reaction zones exhibits a linear decrease tendency along the axial height of the 
riser, and the average velocity is presented by the mean value of the velocity at 
the inlet (ui) and the velocity at the outlet (u 2 ), which is calculated as 
described hereinbelow: 



4). Detailed Deduction Procedure : 

Based on the above-mentioned prerequisites or assumptions for 
deduction and the data disclosed in EP0171460 (the height of the prelift zone 
is 4.9 m, the height of the first reaction zone is 14.7 m, the height of the 
conjunct section between the first and second reaction zone is 4.9 m, the 
height of the second reaction zone is 24.5 m, the total height of the riser is 49 
m and the average velocity in the whole riser is 24.5 m/s), the residence time 
in the whole riser is calculated as described hereinbelow: 

4.9m 14.7m 4.9m 24.5m 49m 

— + u i- i u o~ + K±^I + ^r Z= 24^7"s (4) 



It can be seen from the data provided by the Examiner that Uo=9.8 m/s, 
and thus Equation (4) may be written as: 

4.9m 14.7m 4.9m 24.5m 

9T8m7^ + ^ir-9-8m/s + uT+H^^T = 2s ( 5 ) 
InUi-lng.Sm/s 2 



According to the law of conservation of mass, the mass flowrates at all 
cross-sections of the riser are equal, and, in terms of Figure 1, the mass 
flowrate at the inlet of the conjunct section between the first and second 
reaction zones (shown by the subscript "1") and the mass flowrate at the outlet 
of the conjunct section between the first and second reaction zones (shown by 
the subscript "2") are equal, as shown by the following equation: 



p 2 X V 2 = Pl X V 1 (6) 



The general knowledge of the principles of catalytic cracking teaches us 
that the reactions occurring at cross-section "1" and cross-section "2" are 
basically conversion reactions and cracking does not occur in a large scale. 
Accordingly, it may be concluded that the number of molecules at said two 
cross-sections remains unchanged. Since the pressures and temperatures at 
said two cross-sections do not change much, it may be concluded, on the basis 
of the density definition, i.e. equation (2), that the densities at said two 
cross-sections are equal: 

P 2 = Pi (7) 

Equation (6) can thus be written as: 

V 2 = V t (8) 

Upon substitution of the mass flowrate definition, i.e. Equation (3), into 
Equation (8), a new equation is obtained, which is as follows: 

U 2 XTtx(^) = Ul XTTx(^) (9) 

Equation (9) is transformed to obtain the following equation: 




The following monadic nonlinear equation about Ui is thus obtained by 
substituting Equation (10) into Equation (5): 



4.9m 24.5m 



The Examiner believes that the ratio of the diameter of the second reaction 
zone to the diameter of the first reaction zone is 3, which means that x=3 in 
Equation (1 1). The following result is then obtained by solving Equation (11): 



Uj = 183.3m/s (13) 



Obviously, 

u x 183.3m/s 
u 0 9.8m/s 



18.7 (14) 



The ratio of 18.7 is far greater than the generally recognizable ratio of 
less than 3 in the art. 

Even when the ratio of the diameter of the second office zone to the first 
office action is 1.5:1, i.e. x=l.5, The following result is then obtained by 
solving Equation (11): 

«, = 62.8 m/s 

Obviously 

u, 62.8 r A 

- 1 = = 6.4 

u 0 9.8 

This value is even greater than 3 . 

Thus, the Examiner's opinion that "DR 2 : DRi is equal to approximately 3: 
1" is incorrect and impossible. 

5) To further demonstrate that U]/u 0 <3, please see the following deduction 
based on prior art. 

Deduction 1 : 

According to the law of conservation of mass, the mass flowrate at the inlet 
of the first reaction zone (corresponding to cross section "0" in Figure 1) and 
the mass flowrate at the outlet of the first reaction zone (corresponding to cross 
section " 1 " in Figure 1) are equal, i.e. 

P : xV a =P o xV 0 (15) 



The definition of density (equation (2) ) and the definition of the volume flow 
(equation (3) ) are substituted into equation (15), and it produce: 
Pi xMj P 0 X M 0 



Transform equation (16) and produce: 



HI)*©*® a» 



According to the general knowledge of the principles of catalytic cracking, one 
skilled in the art knows that the range of variation of pressure at the inlet and 
outlet of the first reaction zone is narrower than that of temperature, i.e. in 
equation (17): 



Considering equation (18) in equation (17), it is easy to see that: 



Since reactions basically including cracking occur at the bottom of the riser 
(the so-called first reaction zone) after large-molecular feedstock oil contacts 
with a hot regenerated catalyst, the conversion of the feedstock in said zone 
increases rapidly, and a large amount of small-molecular products are 
produced, which is the manifestation of molecular expansion. It can be known 
that according to the typical distribution of catalytic cracking products, the 
ratio of the molar mass at the inlet and the outlet of the first reaction zone in 
the equation (19) is usually less than 3. The patent Examiner can find out that 
the distribution of many catalytic cracking products produces this ratio by 
calculation. Hereinbelow is a specific example for the Examiner's reference in 
the calculation of this ratio. 




(18) 



The calculation process of the gas composition and the molar mass thereof 
at the inlet of the first reaction zone is: 





Feedstock Oil Gas 


Atomized Vapour 


Mass Composition, wt% 


90 


10 


Molar Mass, g/mol 


500 


18 


Mol 


90g 


10g 


SOOg/mol 


18g/mol 


Total Mol 


90g lOg 
500g/mol 18g/raol 


Molar Mass of the Mixed Gas, 
g/mol 


100g 

90g lOg " 135 - 95 s/ m ol 
500g/mol 1 18g/raol 



The calculation process of the gas composition and the molar mass thereof 
at the outlet of the first reaction zone is: 





Dry Gas 


Liquefied 
Gas 


Gasoline 


Diesel 
Fuel 


Cycle 
Stock and 
Slurry Oil 


Atomized Vapour 


Mass 

Composition, 
w% 


2.7 


18 


36 


18 


15.3 


10 


Molar Mass, 
g/mol 


20 


50 


100 


200 


500 


18 


Mol 


2.7g 


18g 


36g 


18g 


iS.3g 


10g 


20g/mol 


SO g/mol 


lOOg/mol 


200g/raol 


500g/:raoJ 


18g/mo] 


Total Mol 


2.7g 18g 36g 18g lS.-3g lOg 
20g/raol 50g/mol lOOg/mol 200g/mol. 500g/mo] ' 18 g/mol 


Molar Mass of 
the Mixed Gas, 
g/mol 


100g 


2.7g , lSg , 36g , 18g , 15.3g , lOg ' - 65..31g/mol 


20g/moi 1 SOg/mol ' lOOg/mol ' 200,g/mol ' 50Og/mol ' i8g/mol 



Thus, the ratio of the molar mass of the mixed gas at the inlet of the first 
reaction zone and at the outlet of the first reaction zone is: 

M 0 _ 135.95g/mol 
Mj 6S.31g/mol 



i.e. in equation (19) 



Deduction 2: 



I have provided three known references in my previous Declaration. They 
also prove my above viewpoint. 

1) Figure 2 which illustrates the distribution of conversion along the axial 
direction of the height of the riser (In-Su Han, Chang-Bock Chung, Chemical 
Engineering Science 56, 2001, page 1954, see Annex I in my previous 
Declaration). 



1 st reaction zone 2 nd reaction zone 




Riser height (m) 

Fig 2: Distribution of conversion along the axial direction of the height of the riser 

It is generally accepted in the art that the ratio of the velocity at the outlet 
of an iso-diameter riser( i.e. the so-called first reaction zone) to the velocity at 
the inlet of the iso-diameter riser( i.e. the so-called first reaction zone) is less 
than 3. This is exemplified by the following two examples: 

First Example : Figure 3 illustrates the distribution of velocity along the 
axial direction of the height of the riser (In-Su Han, Chang-Bock Chung, 
Chemical Engineering Science 56, 2001, page 1978, see Annex II in my 
previous Declaration). It can be seen that the ratio of the velocity at the outlet of 
the first reaction zone to the velocity at the inlet of the first reaction zone (i.e. 
Uj/uq in Figure 2) is less than 2. 




Fig. 3 Distribution of velocity along the height of the riser 

Second Example : It can be readily seen from the reference Process 
Calculation and Technical Analysis of Catalytic Cracking (written by Cao 
Hanchang, Hao Xiren and Zhang Han, see Annex III, and the translation 
thereof, Annex IV in my previous Declaration), for example, the calculation 
results on page 242 (see also pages 5 and 6 of Annex IV) show a velocity of 
11.45 m/s at the outlet of the riser, a velocity of 5.78 m/s at the inlet of the 
riser and a ratio of the former to the latter of 1 .98. 

9. Conclusions: 

(1). The drawing in the reference (EP0171460, Kmecak et al) cited by the 
Examiner is a schematic drawing. An actual size can not be derived form this 
drawing. 

(2). The velocity inside the riser of 18~31m/s or the residence time of 2s, 
mentioned by Kmecak et al , is much contradictory to the Examiner's opinion 
that the ratio of the diameter of the second reaction zone to that of the first 
reaction zone is 3, even 1.5. 1 respectfully deem that the Examiner's comments 
are incorrect. 

10. I further declare that all statements made herein on personal knowledge are 
true and that all statements made on information and belief are believed to be 
true; and that these statements are made with the knowledge that willfully or 
intentionally false statements made by declaration to the United States Patent 
and Trademark Office are punishable by fine or imprisonment, or both, under 1 8 



U.S.C. §1001; and that such willfully or intentionally false statements may 
jeopardize the validity and/or enforceability of the patent application with Appl. 
No. 09/553,990 and any United States Patent issuing therefrom. 

Signed at /^j''^ this I ^ day of M^oj ^ 2 008. 

Dr. Xu Youhao 
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Abstract 

The aim of this study is to obtain a model that can simulate the performance of an industrial fluid catalytic cracking (FCC) unit in steady 
and dynamic state, and which will subsequently be used in studies of control and real time optimisation. Tn this paper, a dynamic model for a 
R2R type FCC unit is presented. The model includes the riser, the stripper/disengager, the regeneration system and the catalyst transport lines. 
Mass, energy and pressure balances are performed for each of these sections. 

Simulation results for steady state are presented and compared qualitatively to those obtained from previous FCC models. The dynamic 
behaviour of the system is explored through two perturbations in open loop, one on the fresh feed flow rate and one on the air flow rate to 
the first regenerator. The results illustrate the consistency of the model and are in agreement with what has been observed in studies available 
in the open literature. 
€■ 2006 Elsevier Ltd. All rights reserved. 

Keywords: Fluid catalytic cracking; Mathematical modelling; Dynamic simulation; Nonlinear dynamics 



1. Introduction 

Fluid catalytic cracking (FCC) is one of the most important 
processes in a refinery, and is therefore sometimes referred to 
as the heart of the refinery. Several studies have been made 
on the modelling, simulation, kinetics, multiplicity of steady 
states, chaotic behaviour, on-line optimisation and control of 
FCC units. However, there are still large areas to be examined 
due to the complexity and to the economical importance of this 
process. 

The initial attempts to model FCC units where primarily 
oriented towards the kinetics of cracking. The first and most 
widely used kinetic models are the three lump model proposed 
by Weekman and Nace (1970), which mainly focused on feed- 
stock conversion and gasoline selectivity, and a 10 lump model 
that introduced a more chemical description of the heavier 
fractions, also proposed by Weekman and co-workers a few 
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years later (Jacob et al., 1976). An extended kinetic model that 
uses 18 lumps to distinguish several chemical groups within 
each hydrocarbon fraction has been published by Pitault et al. 
(1994). 

More detailed kinetic modelling approaches have also 
been applied to catalytic cracking. The "Structure Oriented 
Lumping" approach, described by Quann and Jaffe (1992, 
1996) and Christensen et al. (1999), uses an extensive net- 
work of reactions between lumps that are represented as an 
assembly of molecular building blocks. In more recent studies, 
detailed kinetic models for the catalytic cracking of industrial 
vacuum gas oil (VGO) feeds have been developed based on 
the single events approach (Moustafa and Froment, 2003). 
This modelling technique aims at retaining the full detail of 
the reaction pathways of all individual feed components and 
reaction intermediates by describing the reaction network in 
terms of elementary steps. It therefore requires a very detailed 
characterisation of feeds and products, which is generally not 
available in industry, and is computationally more demanding 
than lumped models. Semi-lumped models that consider the 
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reaction types occurring in catalytic cracking have also been 
proposed, in which the products and reactants are lumped by 
characteristics such as the carbon number (Pinheiro et al., 1999; 
Carabineiro et al., 2004). Unfortunately, studies concerning this 
type of models have been made mainly for model molecules and 
like the single event models they require good characterisation 
of the reactants and products. 

In FCC modelling, studies of the regenerator hydrodynamics 
and steady state behaviour have also been conducted early on 
(Errazu et al., 1979; De Lasa and Grace, 1979). The combustion 
phenomena in the regenerator, such as the CO post-combustion 
reaction (Morley and De Lasa, 1988) and the prediction of the 
CO2/CO ratio in the flue gas (Weisz, 1966), have been a matter 
of interest for investigators. 

Several authors have presented steady and unsteady state 
models that describe the riser reactor-regenerator system of a 
FCC unit. Among these are works by McFarlane et al. (1993), 
Arbel et al. (1995a), Han and Chung (2001a), Elnashaie et al. 
(2004) and Secclii et al. (2001), which will be described in 
more detail. 

McFarlane et al. (1993) published a well-detailed model 
based on the Exxon model IV that covers most parts of a Model 
IV type FCC unit. However, the model uses an oversimplified 
cracking kinetic model and it also lacks a description of com- 
bustion reactions. Arbel et al. (1995a) developed a model to 
predict both steady and unsteady states of the current genera- 
tion of FCC units. For the riser reactor, their model assumed 
pseudo steady state, plug flow, no slip velocity between the gas 
and catalyst particles, constant superficial velocity and adia- 
batic operation. The cracking kinetics were described by means 
of the 10 lump model proposed by Jacob et al. (1976). The 
regenerator was modelled according to the two-phase model 
described in Kunii and Levenspiel (1990). Arbel et al. (1995b) 
used this model in their study on the multiplicity of steady 
states and in control applications. 

Elnashaie et al. (2004) have recently presented a dynamic 
model of the riser reactor-regenerator system that is similar to 
the one developed by Arbel et al. (1995a). In their paper, they 
mainly focus on bifurcation and multiplicity of steady states in 
the FCC unit that arise from the highly exothermic combustion 
reactions and back mixing in the regenerator. 

Secchi et al. (2001) presented a model for a stacked UOP 
FCC unit where the riser is described as a plug flow reactor 
in dynamic state, however, with no overall mass accumulation. 
Slip velocity between gas and catalyst was not considered. Two 
hydrodynamic models were proposed for the regenerator: a 
CSTR and a two-phase theory model. After comparison of these 
two approaches, the authors concluded that the CSTR is not able 
to correctly predict the CO/CO2 ratio. On the other hand, Errazu 
et al. (1979) and Vale (2002) also compared steady state results 
with a two-phase theory model for the regenerator to those 
obtained with a one-phase model that represents the regenerator 
as a CSTR. They both concluded that there were no significant 
differences between the predictions of the two models, which 
is in disagreement with the conclusions of Secchi et al. (2001). 

A quite comprehensive model of the FCC unit was presented 
by Han and Chung (2001a). Their model describes a FCC unit 



comprising a riser, stripper/disengager, regenerator and catalyst 
transport lines with slide valves. Besides the mass and energy 
balances for the riser, the authors also wrote the momentum 
equations for the gas and the catalyst particles, which allows 
determining the velocity profiles along the riser reactor and ac- 
counts for the slip velocity between gas and catalyst particles. 
The authors claim that a deviation in gas oil conversion as high 
as 9.6% can be obtained by including and excluding the mo- 
mentum equations in the riser model (Han and Chung, 2001a). 
In their second paper (Han and Chung, 2001b), however, they 
present velocity profiles with a slip velocity below 0.25 m/s, 
corresponding to a slip factor close to 1. This is clearly in dis- 
agreement with typical slip factors in the riser of about 2, which 
have been used by many other authors (Das et al, 2003; Fligner 
et al., 1994; Malay et al., 1999). 

The existence of high slip factors is attributed by Fligner et 
al. (1994) to the formation of clusters along the riser reactor. 
A cluster is an aggregate of particles moving together with the 
same velocity and has a higher "free falling velocity" than 
the terminal velocity of individual particles. Hence, when the 
amount of clusters increases, the slip factor between gas and 
catalyst will also increase. According to Parssinen (2002), the 
reported size of FCC catalyst clusters is between 2 and 15 mm. 
In the dilute region of the riser, however, the probability of 
cluster formation is smaller due to the higher velocities and the 
slip factor decreases significantly in this region. 

Currently, several types of FCC units are in operation, using 
various designs. Hence, not all models available in literature 
can easily be applied for every type of unit. In this study, a dy- 
namic model is presented for an R2R unit comprising a riser, 
a stripper, a disengager, two standpipes and a regeneration 
system composed of two regenerators connected by a lift. An 
effort was made to develop a model that is as simple as possi- 
ble at the computational level without disregarding important 
characteristics of the FCC unit, such as its hydrodynamics and 
a cracking kinetic scheme able to predict the major cracking 
products: decanted oil (DO) (unconverted vacuum gas oil feed); 
light cycle oil (LCO), gasoline (GLN), liquefied petroleum gas 
(LPG), fuel gas (FG) and coke (CK). The riser hydrodynamics 
accounts for the molar expansion for the gas and considers a 
force balance on the catalyst particles, which results in a slip 
velocity between catalyst and gas. For the regenerators, a sim- 
ple CSTR hydrodynamic model was used; a previous study on 
R2R units (Vale, 2002) showed the good performance of this 
model, since there were no significant differences between the 
steady state behaviour of a two-phase theory model and a CSTR 
model. 

The FCC model was implemented in FORTRAN using a 
modular approach. A qualitative analysis of steady state results 
and of the dynamic behaviour has been made, showing promis- 
ing results for future use of this model in studies of advanced 
control and real time optimisation. 

2. Process description 

The R2R resid FCC process, initially devised by Total 
(Dean et al., 1982a,b) and now licensed by Axens/IFP and 
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Stone & Webster (Marcilly and Bonifay, 1996; Letzsch and 
Jackson, 2004), consists of a riser reactor, a stripper and two 
regenerators. The two-regenerator R2R technology allows to 
process feeds with a high residue content, as more coke can be 
burnt. In this configuration, the first regenerator (regenerator 
1) acts as a mild pre-combustion zone that removes 40-70% of 
the coke on the catalyst. The partially regenerated catalyst with 
less than 0.5 wt% coke is then air-lifted to the elevated second 
regenerator (regenerator 2), where complete regeneration is 
achieved with slight air excess and under a low steam partial 
pressure in order to minimise catalyst deactivation (Gauthier 
et al., 2000). A schematic of a typical R2R resid FCC unit is 
shown in Fig. 1 . 

The regenerated catalyst flows through the standpipe from 
the regenerator to the bottom of the riser. After a short re- 
acceleration zone to stabilise the catalyst flow, preheated finely 
atomised oil feed is injected onto the hot catalyst. Upon con- 
tacting the catalyst, the feed vaporises and cracking starts. 
Vaporisation and cracking generates a significant expansion 
that lifts the catalyst up the riser. As the cracking reactions 
proceed throughout the riser, heavy hydrocarbons are cracked 
into desired products, but large amounts of coke, 4-8 wt% of 
the feed, are simultaneously deposited on the catalyst, reducing 
its activity. 

At the end of the riser, a riser termination device rapidly 
separates hydrocarbon vapours and catalyst particles to reduce 
further thermal and catalytic cracking. This section is com- 
posed of a stripper, in which steam is used to remove most 
of the entrained hydrocarbon vapours in a counter-current 
dense phase stripper with multiple steam injections, and of a 
disengagement vessel with cyclones that ensure an efficient 
separation of the entrained catalyst particles. The stripped 
spent catalyst is then introduced on top of the first regener- 
ator fluidised bed, where the hot flue gas provides ultimate 
stripping. 



In the fluidised bed regenerator, the catalyst activity is re- 
stored by combusting the coke in air. Besides restoring the cat- 
alyst activity, the combustion process also raises the catalyst 
temperature, thus providing the necessary energy to the vapor- 
isation of the feed and the cracking reactions in the riser. Tnside 
the regenerator, two distinct regions can be observed: the dense 
region and the dilute region, frequently termed freeboard. The 
dense region concentrates most of the solids and is the phase 
where most of the combustion occurs. At the surface of the 
dense fluidised bed, bubbles erupt, ejecting particles into the di- 
lute phase freeboard, where the solids concentration decreases 
sharply with height towards a very low catalyst concentration 
(De Lasa and Grace, 1979). 

The transport of catalyst particles from the stripper to the 
regenerator and from the regenerator to the riser bottom is 
achieved by means of two standpipes. Smooth operation of 
these standpipes is critical to an FCC unit, since the standpipes 
link riser and regenerator. It is also important to note that the 
catalyst flows to the riser and to the regenerator are controlled 
by a slide valve positioned in each standpipe. The opening of 
these valves controls the catalyst circulation rate through the 
entire FCC unit. 

3. Mathematical model 

The mathematical model presented here describes the steady 
state and dynamic behaviour of a R2R industrial unit. The 
main assumptions are instantaneous and complete feed vapor- 
isation, pseudo-steady state for the riser, which is modelled as 
an adiabatic plug flow reactor with a six-lump kinetic model, 
a CSTR dynamic model for the stripping/disengagement 
vessel and a CSTR dynamic model to describe each regen- 
erator. The lift between the regenerators is also modelled as 
a plug flow reactor in pseudo-steady state. The catalyst cir- 
culation rate through each standpipe is determined by the 
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pressure drop across a slide valve and is given by a valve 
equation. 



3.1. Feed vaporisation 

According to Ali and Rohani (1997), the inlet zone is the 
most complex part of the riser, since there is high turbulence, 
flow inhomogeneity and large temperature and concentration 
gradients. However, the same author also indicates that the feed 
takes only about 3% of the mixture residence time to vapor- 
ise completely, which justifies the assumption of instantaneous 
vaporisation. 

As the vaporisation is assumed to be instantaneous and com- 
plete, the energy balance will be written in the steady state with- 
out any heat transfer resistances. The change in feed enthalpy 
during vaporisation is given by the correlations in Montgomery 
(1 993) that depend on the feed characteristics and can be found 
in the Appendix. 



3.2. Riser 

The riser hydrodynamics is modelled as a plug flow reactor 
with a one-dimensional model that only considers axial vari- 
ations of the variables and diffusional and turbulent transport 
of mass, momentum, and heat are neglected. An hybrid model 
was used to describe gas and solid velocity along the riser. The 
gas phase was modelled as a continuum phase and the continu- 
ity equation that expresses molar expansion was used to obtain 
the gas velocity profile. For the solids phase, it was considered 
that all FCC particles are identical and Newton's second law 
of motion was used to write a force balance over a particle. 

For concentration and temperature; both internal and external 
mass- and heat-transport resistances were neglected, and these 
variables are modelled considering a pseudo-homogeneity as- 
sumption. 

The cracking kinetics is described by a six-lump model for 
gas oil cracking adapted from the model first presented by 
Takatsuka et al. (1987) with a deactivation function depend- 
ing on the catalyst coke content (Vale, 2002). The lumps con- 
sidered are VGO + DO(> 360°C); LCO (220-360 °C); GLN 
(C5-220°C); LPG (C3 and C4); FG (H2, CI, C2 and H2S) 
and coke. The reactions between these lumps are presented 
in Fig. 2. Similarly to other authors (Ancheyta-Juarez et al., 
1999; Corella and Frances, 1991 ; Hagelberg et al., 2002) gaso- 
line and gases conversion to coke were not considered since 
the kinetic constants for these reactions are negligible in com- 
parison to the other reactions. 

Except for the second order VGO cracking reactions, all re- 
actions are first order. The kinetic parameters were determined 
from pilot plant data. 

According to this kinetic model, the rate of the reaction of 
lump i towards lump j is given by 




(adapted from Takatsuka et al., 1987). 

with $ the deactivation function, which is given by 

* = exp(-a7^ al ). (2) 

Mass balances for each lump and for the coke content on the 
catalyst are given by 



(5) 



The coke deposited on the catalyst can be of two different na- 
tures: catalytic coke that results from secondary catalytic reac- 
tions and CCR coke which results from thermal decomposition 
of heavy compounds that are present in the feed. According to 
Arbel et al. (1995a), CCR coke deposits almost immediately 
on the catalyst surface. However, since this coke is not directly 
formed at a catalytic site and is mainly deposited on the cata- 
lyst matrix, it is considered that it does not contribute to cata- 
lyst deactivation. From the total amount of CCR present in the 
feed, only a fraction xqcr is converted into coke. The follow- 
ing balance then gives the total amount of coke on the catalyst 
in the riser: 

-fck, rs --- Y°g RS + Y-f* (6) 
with 

F9 CR Xrrr> 

(7) 



Gas molar expansion as well as slip velocity between gas and 
solids is accounted for in the riser model. Formation of clusters 
of catalyst particles has also been considered, since several 
papers show evidence of this phenomenon for FCC catalysts 
(Fligner et al., 1 994; Parssinen, 2002). The size of a cluster was 
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determined by using the correlation from Xu and Kato (1999) 
that can be found in Appendix. The molecular weights of the 
lumps needed to predict gas molar expansion were taken from 
Pitault et al. (1994). 

To derive the velocity profile for the gas phase the continuity 
equation as a molar balance was written as already mentioned. 



" s at + dz ~ {l £s) ^M,„r 



(8) 



tog 
dz : 



C, V ' C, dz z s dz ' g dt ' 



For the gas phase it is assumed that hydrocarbon vapours be- 
have like an ideal gas and the ideal gas law is used for total 
concentration calculations. 



(12) 
(13) 



6Q 1 QP RS 
dz ~ RT dz 



T dz 



For the solid phase, a force balance was done over a cluster of 
FCC particles. The movement of a FCC cluster is then given by 
the equation of Newton (Pugsley and Berruti, 1996; Sabbaghan 
et al., 2004): 



ndl. dv r pAv s -v c f ndr. ndl 



The drag coefficient, Cd, for the gas-solid interphase friction 
is given in Appendix. 

Considering the equality: 
dv c _ dz dv c . dv c _ Qv c 

~oT ~ l! * o7 ^ Hi 

Eq. (14) can be rearrange 
the riser for (he solid phas 



o give, the velocity profile along 



= C D ~ 



(16) 



For the energy balance the riser was considered as an adiabatic 
phig flow reactor without heat transfer resistances between gas 
and solid phases. An overall cracking enthalpy expressed per 
unit mass of VGO converted covers the overall heat consumed 
by the endothermic cracking reactions. Due to the assumption 
of pseudo-homogeneity for the energy balance both gas and 
solids share the same temperature (Das et al., 2003). 



87bs , P B C P ,8+F C C P , C 87» s _ 



■■e.p t C„. + z cPc C„, 



SDQ, 



(17) 
(18) 



Pressure profile along the riser is given by the total momen- 
tum balance for both gas and solid phase (Arastoopour and 
Gidaspow, 1979). 

_ = _ (£cPc + esPg)g ... —[pgtgvj; 

+ p c {\ - s s )vl] - F w . (19) 

Pugsley and Berruti (1 996) showed that the total pressure drop 
per unit length along the riser in fully developed flow can be 
approximated by the pressure drop caused by hydrostatic head 
of solids and that the acceleration effects are important only at 
the base of the riser. Frictional effects due to solid friction and 
gas friction with the riser walls, F w , can then be neglected. 

9P»s 



(20) 



The void fractions of the si 
following definitions: 



" v c p c Qn S ' 



:: 1 - S, 



d and gas phases are given by die 

(21) 
(22) 



For dynamic simulation purposes the riser is considered to be 
in pseudo-steady state, since it has a residence time of the 
order of a few seconds (2-5 s), while the regenerators and the 
stripper have characteristic times in the order of a few minutes. 
It was then considered that the regenerators and stripper vessels 
control the transient response of the system reactor-regenerator 
when perturbations to the system occur. 

3.3. Stripper/ disengager 

In recent years many developments have been made in the 
design of riser termination devices in order lo achieve a bet- 
ter separation of the gas and solid phases and to reduce over- 
cracking of gasoline into lighter components. In this work, it 
will be assumed that the separation of gases and solids is fast 
and efficient so there are no reactions occurring in the disen- 
gager/stripper section, which will be modelled as a continuous 
stirred tank (CST) in unsteady state. The properties of (he cat- 
alyst bed will be determined by considering that the fluidised 
stripper bed is at incipient fluidisation. Since the catalyst level 
in the stripper is an important variable for the control of the 
FCC unit, the variations of the stripper bed level will be ac- 
counted for (Eq. (24)-(26)). 

It is also assumed that a percentage of hydrocarbons remains 
adsorbed/occluded in the catalyst pores after stripping. This is 
the so-called cat-to-oil coke (Montgomery, 1993). An empirical 
stripping function based on pilot plant data, given in Appendix 
is proposed for predicting cat-to-oil coke. 

Coke mass balance: 



or ::=: Wc.sv 



(23) 
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3Wr,ST — Pc e cQsTVLsT, 
HLst _ f'c.RSout - ^t.STout 



Gas inventoiy: 



^S,RSout + F s - J^.STout- 



(24) 
(25) 
(26) 

(27) 



F« s 



(28) 



The catalyst flow rate leaving the stripper towards the regener- 
ator is controlled by a slide valve in the spent catalyst standpipe 
and the pressure difference between these two vessels: 



APsp s - /"sTb^cn -4- p c s c Ls? s sin d& s g - J^gi 



(29) 
(30) 

'Hie pressure in the disengagement vessel is calculated by as- 
suming ideal gas behaviour for the hydrocarbon vapours: 



V g ,ST = VST 



W c ,sr 



(31) 
(32) 



The pressure at the bottom of the disengager/stripper section 
is higher tban at the top, due to the static head exerted by the 
catalyst. This is given by 



The energy balance is written by assuming thermal equilibrium 
between the outlet streams, neglecting any heat transfer with 
the exterior. The heat of desorption will also be neglected. 

fl(W c ,STCV t TsT + W g ,STC p /isT) 

di 

= ^s.RSontHg.RSout 

+ Fv,StH.,,ST-t- Jt-.RSouiH^RSont 
~~ •fg.STovil.Hg.SToiU — ^c.SToiuHr^STouti 

H 4 = C p , x (71 - T„f) 
3.4. Regenerators 



(34) 
(35) 



FCC regenerators consist of large fluidised bed r 
where the strongly exothermic combustion reactions take place. 
In this study, the freeboard will be neglected and the dense 
phase will be modelled as a CSTR. 



Coke is considered to be composed of carbon and hydro- 
gen, although in practice sulphur and nitrogen are also present 
in small quantities. Moreover, no distinction is made between 
different types of coke with respect to its combustion. This 
last assumption might be subject to some discussion, as we 
consider that gaseous hydrocarbons entrained with the catalyst 
(soft coke) burn in the same way as hard coke, which is a solid 
deposited on the catalyst. 

The combustion kinetics considered in this work was taken 
from previous studies at EFP (Daly et al., 2001; Vale, 2002; 
Zwinkels and Nougier, 1 997). The combustion of carbon can 
be written as the following overall reaction: 



The gas flow rate at the disengager outlet is determined by the 
opening of the valve to the main tractionator, which depends 
on the pressure in the stripper and in the main fractionator. The 
gas flow rate is therefore given by the following equation: 



C + 



CO,+ 



I 



CO. (36) 



The combustion of carbon produces both CO2 and CO in a 
proportion of r>. The CO2/CO ratio a is very important, since 
it affects both the mass and energy balances in the regenerator. 
Indeed, the heat of combustion to CO2 is almost 3 times the 
heat of combustion to CO, which means that the CO2/CO ratio 
controls the heat balance in the regenerator (Arbel et al., 1 995a). 

Several correlations have been derived to calculate a, but 
the one proposed by Arthur (1951) is the most quoted in the 
literature. Arthur's correlation states that a decreases exponen- 
tially with temperature. Experimental work done at IFP (Daly 
et al., 2001; Zwinkels and Nougier, 1997), however, has shown 
that the carbon combustion rate and the a ratio are affected by 
other variables than the temperature, such as water concentra- 
tion, coke conversion and metals content. 



c + ±o 2 — * CO, 



rc = (V Al y + V W )Y C , 

X p of + k C2 ex P (-^ t /i?7)P H2 o ' 



k™ l em~E™ l /RT)P™ 

d0.11 

'jh 

= *£exp(-Ej£//?7)/(Xc)' 



(37) 
(38) 

(39) 
(40) 
(41) 
(42) 

(43) 
(44) 
(45) 
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The subscripts dry and wet represent the partial contributions 
due to dry and wet conditions during coke combustion, while 
/ (Xc) is a function of carbon conversion. The kinetic constants 
were determined at IFP (Daly et al., 2001) and are given in 
Table 4. 

In most models, the hydrogen combustion reaction is gener- 
ally disregarded. However, according to Faltsi-Saravelou and 
Vasalos (1991), this reaction has to be considered since it has 
a significant thermal effect. In this study, the rate expression 
for the hydrogen combustion reaction obtained by Wang et al. 
(1986) will be used. 



(47) 
(48) 



The carbon combustion reactions produce simultaneously CO; 
and CO. However, the CO produced undergoes further oxida- 
tion to CO?, through the so-called after-burning reactions. The 
oxidation of CO to COt can be of two different natures: het- 
erogeneous (catalytic) or homogeneous combustion. 



CO + 



"— > co 2> 



(49) 



co + ±02-^ ah. (50) 

Tn this work, the following rate expression for the heteroge- 
neous combustion, determined at IFP by Daly et al. (2001), will 
be used: 



p o.4y p -( 



The homogeneous conversion of CO to CO? is a free radical 
reaction which is inhibited by the presence of surfaces such as 
the catalyst (Arbel et al., 1995a). The rate expression used in 
this study was obtained by Zwinkels and Nougier (1997). 

8 k*ap(.-E* al /RT) ofi u5 
r i :::: ~ fzi >'o, >'co v H ,o- 

The CSTR regenerator model can now be described as follows: 
Gaseous species molar balance: 

. 0C <.RO _ v « ... V J 

\,RGirj 'Vuil.K 



(52) 



Vbed- 



lt X C/.RG 



Carbon and hydrogen mass balance: 

Wc.RG B7/,R G _ FcSLGm y . ^;,ROout y 

M,,,,/ <H M w j ,,1RGm M Wti '' RG 

+ V bedF . cP , EC'1 ' :: = C or H. (54) 

j 

Catalyst inventory: 

OLrg JV.RGm ~ -Fc.RGout 



Gas inventory: 

= Fg.RQin - F ff ,RGont 

+ ^(Fc^nYi^-^iFa^Yi^G). (56) 

Energy balance: 
8(W,C ; , ! Trg + W.,r G C /j? 7 rg ) 



(57) 



0/ 

= ^ s ,RGmHg,RGm + ^c.RGin Hc.RGin + 0,° 
- ^ s ,RGoutH s ,RGoitt - /v,RGoutHc,RGont- 

Pressure: The pressure in the dilute phase of the regenerators 
is given by the ideal gas law, while for the pressure at the 
bottom, the hydrostatic pressure exerted by the catalyst bed at 
the bottom of the regenerators is considered. 

flMW, :::: p s,*G + p c tc L&Qg. (58) 
The catalyst flow rate that passes through the standpipe that 
links the second regenerator to the riser is then given by 

F,.RG2out = F,:,SP R = * U ,SP„ (x v )/AP S P~, (59) 

APsp* = Prgz + P c £cLsp s sin 6 S p K g - PRs b(lltom . (60) 

The gas flow rate exiting from the regenerator depends on the 
pressure difference between the flue gases coming out from the 
regenerators and the pressure inside the regenerator. 



f^RGout = fc«,RG (*.)) V^.RG - 



(51) 3.5. Lift 



The two regenerators are linked by a lift that transports cat- 
alyst from the lower regenerator to the upper one. Since the 
residence time of the gas and catalyst in the lift is of the order 
of a few seconds, pseudo-steady state will be considered as for 
the riser. It was assumed that both gas and catalyst are in plug 
flow with combustion reactions occurring. 

Gaseous species molar balance: 



Carbon and hydrogen molar balance: 



•'W>,.A..ifi W) . ™i 

Energy balance: 
B7Lift 



,, „ "'Lift . „ „ O/Lift 

= "flLift {*g Ytf AH*) + &cPe E(rjAHpj . 
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Since the molar expansion is negligible in the lift, no velocity 
profiles were calculated and the porosity in the lift is given by 
an empirical correlation available at IFF (Vale, 2002). 



2 1 06C c .exp(~0.251tf i: ) 



where G c is catalyst flux and u 



(65) 

is the gas superficial velocity: 





32.11 




Stripper 


8.0 


3.6 


Disengager 


8.5 


5.1 


Regenerator I 


1 3.7 


6.4 


Regenerator 2 


11.8 


5.9 


lift 


20.64 


0.8 


Spent catalyst standpipe 


8.5 


0.8 


Regenerated catalyst standpipe 


7.5 


0.8 



4. Numerical and computational methods 

The mathematical model of an industrial FCC unit presented 
in this paper corresponds to a large system of differential and 
algebraic equations. To solve this problem, the model was im- 
plemented in FORTRAN using a modular approach. 

Each module corresponds to a section of the unit that is 
solved at each time integration step by an appropriate numerical 
method. The solvers used in this study were: 

• ZERO IN — used to find the roots of nonlinear- algebraic equa- 
tions based on Brent's algorithm (Brent, 1971). 

• LSODE — used to solve systems of ordinary differential equa- 
tions (ODE) based on the Gear's algorithm (Hindmarsh, 
1980; Radhakrishnan and Hindmarsh, 1993). 

» DASPK— used to solve systems of differential algebraic 
equations (DAE) based on the Petzold-Gear's BDF algo- 
rithm (Li and Petzold, 1999). 

Initial and boundary conditions were estimated through steady 
state simulation for a base case of typical industrial operating 
conditions. The operating conditions for this base case are pre- 
sented in Table 1 . The R2R dimensions used for simulation are 
presented in Table 2 and feedstock and catalyst properties are 
listed in Table 3. Finally, the kinetic data used for simulation 
is listed both for cracking reactions and combustion reactions 
in Table 4 . 



Tiible 1 




R2R base ease operating conditions 




Base ca.se operating conditions 


Fresh feed flowrate (kg/s) 


60.82 


Fresh feed temperature (K) 


450.15 


Steam fiovwate entering the riser (kg/s) 


3.58 


Steam temperature (K) 


592.25 


Cata!yst-to-oiJ ratio, C/O (kg/kg) 


6.9 


Stripping steam (kg/s) 


2.08 


Stripping steam temperature (K) 


592.25 


Air to regenerator 1 (kg/s) 


30.71 


Air to regenerator 2 (kg/s) 


11.19 


Air to life (kg/s) 




Air temperature (K) 


428.15 


Pressure in the main i'racr.ionator (Pa) 


I.I x I0 ; ' 


Pressure of the flue gases (Pa) 


1.1 x 10 5 



Table 3 

Feed and catalyst properties 



Feedstock properties 

API (oj 20.0 

Watson characterisation factor, A'-w (dimensionless) 11,9 

Comadson Carbon residue, CCR (wt%) 5.78 

Fraction of CCR converted in coke, xqcr 0.31 

Coke composition CI-Io.8 

Specific heat, C pKC (J/(kgK)) 1255.5 

Heat of reaction, Aff a1t (I/kg VGO) 4.2.5 x 10 s 

Average particle diameter, d„ (in) 7.4 x 10~ 5 

Ouster diameter, da (mi 2.0 x 10 - 1 

Density, p K (kg/m J ) 1450 

Specific heat, C pc (I/(kg K)) 1 1 97.5 



5. Simulation results and discussion 

In this section, die simulation results for the steady and 
dynamic states are presented and discussed. Fig. 3 shows the 
influence of the steady state catalyst-to-oil ratio (COR) on sev- 
eral process variables. It can be seen that increasing the COR 
leads to an increase in riser outlet temperature. Since a higher 
flow rate of catalyst enters the riser for a given hydrocarbon 
flow rate, this also contributes to a higher conversion of the 
feed as the curve for VGO yield in Fig. 3b shows. Despite a 
slight increase in coke yield (CK in Fig. 3b), the increase in 
COR leads to a lower coke content on the spent catalyst. This 
decrease in catalyst coke content leads to lower regenerator 
temperatures (RGi and RG2, Fig. 3a), in spite of a slight in- 
crease in overall heat release in the regenerators. Moreover, 
the residence time of the catalyst in the regenerators is lower 
too, due to the higher catalyst circulation rate. These temper- 
ature evolutions with COR are in agreement with the results 
presented by Malay et al. (1999). However, no maximum was 
found in the regenerators' temperature profiles, as presented 
in the works of Arbel et al. (1995b) and Han and Chung 
(2001b). Nevertheless, it is important to note that in their works 
the regeneration system is different from the one presented 

The effect of COR on VGO conversion and product yields 
presented in this paper is also in agreement with the results 
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is used in this study 



Pre-exponential fa 



Activation energy (J/rnol) 



Cracking kinetic pt 
VGO -> LCO 
VGO -> GLN 
VGO -> LPG 
VGO -> FG 
VGO ->■ CK 
LCO ->• GLN 
LCO — CK 
GLN ->• LPG 
GUN -> FG 
LPG -> FG 



6.59 x 10 6 (in 3 /s/kg) 
1.02xl0 7 (mV-Vkg) 
2.20 x 10 s (nv 1 /5/kg) 

7.35 x 10 s (m 3 /»Ag) 

6.36 x I0 5 fm 3 /s/kg) 
8.44 x 10 5 (s-' )i 
7.81 x 10 5 fs" 1 ) J 
6.95 x 10' (s- 1 ) j 
2.24 xl CP (a- 1 )/ 

i .87 x I0 n (s _1 ) 



146 x 10 3 
193 x I0 3 



Heterogeneous CO oxidatio 
Homogeneous CO oxidatior 



9.45 x 10 2 exp(~1.43 x Vfi/KT) (R-'Pa" 0 - 78 ) 

2.35 x 10 4 exp(-7.92 x itf/XT) (Pa"") 
7.23 x 10 l exp(~1.59 x itf/RT) (s-'Pa"') 
5.I0 x I0 5 exp(-l.77 x itf/RT) (Pa"" 016 ) 
1.21 x 10 3 exp(-2.71 x lO 4 //?'/') (s-'Pa" 01 ) 
2.0Xl - 2.5Xl +X C + 0.1 

4.03 x JO 6 exp(-l .58 x I0 s //?T) (mol/s/k.g cat) 
3.89 x Ur 6 exp(-1.15 x KT*/R7-) (mol/s/kg cat) 
8.45 x 10- 2 expf-2.53 x ]Q S /RT) (mol/s/m 3 ) 




Catalyst-to-Oil (kg/kg) 



Catalyst to Oil (kg/kg) 



Fig. 3. Effect of COR on temperatures and product yields, (a) Regenerator and riser temperature; (b) VGO and product yields. 



presented by Malay et al. (1999). Higher temperatures promote 
higher conversion rates and, at the same time, a higher COR 
also means more catalyst, and hence more active centres avail- 
able for reaction leading again to higher conversions and light 
product yields. 

Fig. 4 shows the profiles of some important process variables 
along the axial coordinate in the riser. It can be seen that most 
of the cracking reactions occur in the initial section of the riser. 
This is an expected result because the highest temperatures and 
lowest catalyst coke content are encountered at the riser inlet. 
The velocity profiles (Fig. 4c) clearly show that there is a slip 
factor between the two phases that decreases along the riser 



for higher velocities. This slip factor between the two phases 
is largely due to the formation of clusters of catalyst particles. 

In Fig. 4c it can also be observed that both gas and solids 
velocities increase along the riser due to a molar expansion of 
about 4 times in the total number of moles in the gas phase. 
This molar expansion makes the solids void fraction decrease 
along the riser, while the gas void fraction increases as expected. 
Fig. 5 shows the dynamic response of the system to a step per- 
turbation in the gas oil feed rate. Decreasing the fresh feed flow 
rate initially results in a steep increase of the riser temperature 
(Fig. 5b, RS) due to a higher COR. This immediately leads to 
a steep increase in conversion (Fig. 5a, VGO). On the other 
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hand, a lower concentration of hydrocarbons leads to a pres- 
sure decrease in the stripper (Fig. 5c, ST), which causes a steep 
decrease in the flow rate of spent catalyst and an increase in the 
regenerated catalyst flow rate (Fig. 5e). The evolution to the 
new steady state is achieved through a pressure balance com- 
pensation caused by a decrease of the catalyst level in regen- 
erator 1 (Fig. 5d, RG1) and an increase of the stripper catalyst 
level (Fig. 5d, ST). Coke-on-catalyst decreases, due to a lower 
feed flow rate (Fig. 5e, coke on spent catalyst) and consequently 
a lower CCR coke contribution. After a small initial increase, 
the temperature in regenerator 1 also decreases (Fig. 5b, RG1), 
since the coke content decreased as well as the residence time 
of catalyst (lower catalyst hold-up). The CO/CO2 ratio shows 
an inverse behaviour of the temperature with steeper decreases 
and increases (Fig. 5f, RGl_CO/C0 2 and RG2_CO/C0 2 ) after 
the perturbation occurred. This is expected since the literature 
shows that a decreases exponentially with temperature (Vale, 
2002). Due to lower temperatures in the regenerators, the tem- 
perature in the riser eventually decreases (Fig. 5b, RS; n ) as well 
as the VGO conversion (Fig. 5a, VGO). The new steady state 
shows near the same VGO conversion due to a higher COR, 
even though the temperatures in the riser and regenerator are 

McFarlane et al. (1993) and Han et al. (2001) also present 
in their papers dynamic simulations in open-loop to feed flow 
rate perturbations, however real comparisons cannot be done 
since all the three works present models for different units. It 



can be seen for example that in our case decreasing the feed 
flow rate leads to lower temperatures and that riser temperature 
shows an initial rise caused by the sudden drop in feed flow 
rate. This is perfectly in agreement with the simulation case 
presented by McFarlane et al. (1993), on the other hand is in 
disagreement with the temperatures transient response for the 
same type of perturbation presented in the work of Han and 
Chung (2001b), where decreasing the feed flow rate leads to 
higher temperatures. 

The second perturbation sets the fresh feed flow rate back 
to its initial value. As can be seen, all variables return to their 
initial steady state. 

A perturbation on the air flow rate to regenerator 1 was also 
made and the dynamic response is presented in Fig. 6. The step 
decrease of the air flow rate to regenerator 1 initially causes 
a steep decrease on the temperatures of the regenerators (Fig. 
6b, RG1 and RG2), since less air leads to lower coke combus- 
tion rates and higher CO/C0 2 ratios (Fig. 6f, RGl_CO/C0 2 ). 
The VGO conversion also presents initially a steep increase 
followed by a steep decrease (Fig. 6a, VGO), since lower coke 
conversions in the regeneration section lead to higher coke 
content on the regenerated catalyst at the riser inlet. 

Besides decreasing the temperature, a lower quantity of air to 
regenerator 1 also leads to a lower pressure (Fig. 6c, RG1). This 
will result in a sudden increase of the catalyst flow rate coming 
from the stripper (Fig. 6e, spent catalyst), since this flow rate 
depends on the pressure difference between the bottom of the 
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Fig. 5. Open loop dynamic response to a step perturbation of -3% in the fresh feed flow rate at 10(K)s, followed by a step perturbation to the initial 
value of fresh feed flow rate al 15000 s. (a) Product yields and VGO conversion; (b) temperature; (c) pressure; (d) catalyst hold-up; (c) catalyst flow and 
cake-on-catalyst; (f) ratio CO/C0 2 and coke conversion. 



stripper and the pressure of regenerator 1 (see Eq. (22)). The 
regenerated catalyst flow rate remains practically unchanged 
(Pig. 6e, regenerated catalyst), since the pressure in the stripper 
(Fig. 6c, ST) and in regenerator 2 (Fig. 6c, RG2) remain almost 
constant. 

To balance the pressure decrease in regenerator 1 the catalyst 
hold-up in the stripper will decrease (Fig. 6d, ST), while the 
bed level in the regenerator 1 will increase (Fig. 6d, RG1). The 
spent catalyst flow rate then returns to a value slightly higher 
than its initial one and this will lead to a new steady state where 



important process variables such as product yields, conversion 
and temperatures are quite similar to the previous one, before 
the perturbation had occurred. 

Increasing the air flow rate to regenerator 1 to its initial 
value will lead the system to the. initial steady state. However, 
the dynamic response of the system is not exactly the inverse 
of the first perturbation in air flow rate at 1000 s, showing 
smoother increases/decreases in some of the process variables 
such as temperatures, coke conversion, CO/COo ratio and VGO 
conversion. 
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Fig. 6. Open loop dynamic response to a step perturbation of -3% in the air to regenerator 1 flow rate at 1000 s, followed by a step perturbation to the initial 
value of air to regenerator 1 flow rate at 15 000s. (a) Product yields and VGO conversion: (b) temperature; (c) pressure: (d) catalyst hold-up; (c) catalyst flow 
and cote-on-catalysi; (f) ratio CO/C0 2 and coke conversion. 



6. Conclusions 

A dynamic model has been developed for the simulation of 
the steady state and dynamic behaviour of a R2R type FCC 
unit. This model includes a riser reactor, a stripper, a disen- 
gager, a two-stage regeneration system and catalyst transport 
lines. 

Studies on FCC units are usually dedicated to traditional FCC 
unit's designs; however residue cracking units are becoming 
more important every day due to the increasing need of pro- 
cessing heavier stocks with an economical lower value. This 



work presents a simulation study on one of this type of units: 
the R2R FCC unit. 

The mathematical model developed allows the prediction 
of the dynamic behaviour of an industrial R2R unit as well 
as of its steady state conditions. The open loop study of the 
dynamics of the R2R unit provides insight into the sensi- 
tivity of the response of tire main variables in the process 
to perturbations to the system. From the simulation results, 
it can be concluded that the model exhibits a behaviour 
that is consistent with the experimental data and literature 
results. 
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An effort was made to assume a compromise between a 
detailed model with heavy computational needs and a simple 
model capable of describing the behaviour of the unit. Future 
work will also include studies of advanced control and real time 
optimisation based on this model. 



Notation 



A Parameter for cluster diameter calculation, di- 

mensioiiless 

Ajj Pre-exponential factor of the reacti on lump i -» 

lump j, s"" 1 orm 3 /(skg) 
d molar concentration of lump i, mol/m 3 

C p specific heat capacity, J/(k:g K) 

Cd drag coefficient for a particle on an infinite 

medium, dimensionless 
Ci mass concentration of lump i, kg/m 3 

CCR Conradson carbon residue in the feed, wt% 

CFFE change of fresh feed enthalpy, Btu/lb 

CK Coke, dimensionless 

d diameter, m 

Eij activation energy of the reaction lump t -> 

lump j , J/mol 

Ec energy parameter of the carbon combustion re- 

action, J/mol 

Ef tj activation energy of the gas reaction j, J/mol 

E a j activation energy of the gas-solid reaction j, 

J/mol 

/ dense-phase fraction, dimensionless 

F mass flow rate, kg/s 

FC Fuel Gas, dimensionless 

F w gas and solids wall friction force, kg/(m 2 s 2 ) 

FF°AVl fresh feed 0 API gravity, dimensionless 
FFE fresh feed enthalpy at riser inlet conditions, 

Btu/lb 

FFK fresh feed Watson characterization factor (K w ), 

dimensionless 

FFT fresh feed temperature or combined feed tem- 

perature, °F 

FFVPE fresh feed vapor phase enmalpy, Btu/lb 
g gravity acceleration, (m/s 2 ) 

C c solidsmass flux, kg/(m ? ' s) 

GLN Gasoline, dimensionless 

Hi enthalpy of phase k in sa:eam i, J/kg 

kc rate coefficient of the carbon combustion reac- 

tion, 

k s . kinetic rate constant of the gas reaction j, 

mols.s _1 .m -3 

kj kinetic rate constant of the gas-solid reaction j, 

mols.s" 1 .kg _i cat 
k v valve flow rating factor, kg/ (s kPa 0-5 ) 

L height, m 

LCO Light Cycle Oil, dimensionless 

LPG Liquefied Petroleum Gas, dimensionless 

M w molar mass, kg/mol 



M2 parameter for cluster diameter calculation, di- 

mensionless 

rijj order of the reaction lump i -> lump j, dimen- 

sionless 

N molar flow rate, mol/s 

P pressure, Pa 

Q n r global reaction heat at standard conditions, J/s 

gi parameter for cluster diameter calculation, di- 

mensionless 

ti rate of formation of lump /, kg/ (m 3 s) 

?ij rate of the reaction lump i -> lump j, 

kg/(m 3 s) 

r? rate of the gas reaction j, mol/(m 3 s) 

rj rate of the gas-solid reaction j, mol/(kg s) 

R universal gas constant, J/(molK) 

Re Reynolds number, dimensionless 

RxT reactor temperature, °F 

t time, s 

T temperature, K 

u superficial velocity, m/s 

v velocity, m/s 

V volume, m 3 

V ( iry partial contribution in dry conditions to carbon 

combustion reaction rate, mol/(kg s) 

V we t partial contribution in wet conditions to carbon 

combustion reaction rate, mol/(kgs) 

VGO vacuum gas oil, dimensionless 

io^? R mass fraction of feed that contributes to cat-to- 

oil coke, wt% 

W inventory, kg 

•*CCR fraction of CCR that converts to coke, 

dimensionless 
x,., valve opening, [0 — 1], dimensionless 

A'c carbon conversion of the coke, kg/kg 

y C k total coke content of catalyst, kg coke/kg 

calalyst 

Y^ catalytic coke content of catalyst, kg catalytic 

coke/kg catalyst 

*dc CR Conradson carbon residue coke content of cata- 

lyst, kg CCR coke/kg catalyst 

Yi component i (in coke) content of catalyst, kg 

i/kg catalyst 

z axial coordinate, m 

Z g gas compressibility factor, dimensionless 

Greek letters 

a, deactivation constant, kg catalyst /kg catalytic 
coke 

y parameter for u>|jj? R calculation, dimensionless 

r volume fraction, m 3 /m 3 

B m f voidage at incipient fluidization, dimensionless 

6sp angle of inclination of the standpipe, 0 

H s dynamic viscosity of the gas, kg/(ms) 

p densityjtg/m 3 

it intrinsic CO2/CO molar ratio, mol/mol 
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stoichiometric coefficient, dimensionless 
heat of cracking per unit mass of VGO con- 
verted, J/kg 

catalyst deactivation function, dimensionless 
cross section, m 2 



Subscripts 




bed 


catalyst bed 


bottom 


refers to the bottom of an equipment 


c 


catalyst 


ck 


coke 


c! 


cluster 


dry 


refers to dry conditions in coke combustion 


8 


gas 




lump or gaseous species 


in 


refers to the entrance of an equipment 


out 


refers to the. exit of an equipment 


P 


catalyst particle 


ref 


reference 




steam 


wet 


refers to wet. conditions it) coke combustion 


MF 


main fractionator 


RS 




ST 


disengager/Stripper section 


RG 


regenerator section 


SP 


standpipes section 


SPa 


refers to regenerated catalyst standpipe 


SPj 


refers to spent catalyst standpipe 


1 


refers to regenerator 1 




refers to regenerator 2 
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Correlations for cluster diameter calculation (Xu and Li, 1998; 
Xu and Kato, 1 999) 



(3333u p g - M 2 )(l ~ tmfKPp ~ P g ) 
(fii - 2M 2 )p p 



(p P - Pg)s r 

Pp L 



(73) 
(74) 
n an infinite 

(75) 
(76) 



Psusp = (' - f)( s g - s cKP r - P g ). 

Correlation for the drag coefficient for a particle o 
medium (Rowe and Henwood, 1961): 

_ f ^(1 + 0.15 x Re p M1 ) Re P < 1000, 
D ~ I 0.44 Re p ^ 1000, 

Rg _ K ~ Vc\ddp g e g 

Correlation for predicting cat-to-oil coke: 

The percentage of feed that contributes to cat-to -oil coke is 
given by the following correlation: 

tugJR = y . C0R) (77) 
7s =e 5 - 2,,3 - ao,44x7 ' OT . (78) 

The total coke concentration on catalyst: in the stripper is then 
given by 

r«*,sr = ^cMS + ?■ (79) 
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Appendix 

Correlations for the enthalpy change of fresh feed during 
vaporisation (Montgomery, 1993): 

FFE === - 13.5 + 0.31FFT 2 + [0.0025(FF°API - 10) 

x(FFT — 100)]+ [0.04(FFK-10)(FFT - 150)] (68) 

FFVPE = 135 + 0.29RxT + 0.000161i?AT 2 
+ [l.3(FF°API-10)] 

+ [0.075(FFK - W)(RxT - 450)] (69) 
CFFE = FFVPE - FFE (70) 
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